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Abstract. A new method of evaluating relative (arm with respect to the interarm disk) star formation rates and 

relative star formation efficiencies, together with spiral arm-amplitudes, as a function of the galactocentric radius, 

using broad-band photometry is derived. The classical method for obtaining star formation rates from Ha photometry 

is discussed, and a new method is derived for diagnosing the possible presence of biased star formation due to different 

initial mass functions in the arms and in the interarm disk. As an example, these methods are applied to the spiral 

Q^' galaxies NGC 4321 and NGC 4254, obtaining their arm amplitudes, relative arm/interarm star formation efficiencies, 

O and relative arm/interarm initial mass functions for each arm, as a function of the galactocentric radius. Both objects 

^ ■ present evidence of massive star formation triggering in the spiral arms consistent with the spiral density- wave theory, 

^ , and a different initial mass function in the arms from that in the rest of the disk, in the sense of favoring a larger 
•I— I , 

fraction of massive stars in the arms. This biased star formation is present in the zones of the arms where there is 
H . . . . . 

^ , triggering of massive star formation, and is then related to, and probably caused, by the density-wave system. However, 
due to this biased star formation, the spiral arms of the galaxies studied do not trigger the formation of a larger total 
mass of stars (of all spectral types) with respect to the interarm disk. 

Key words: galaxies: NGC 4321, NGC 4254 - galaxies: spiral - galaxies: stellar content - galaxies: structure - stars: 
formation 



1. Introduction 

Since the first steps by Lindblad (1959) in the search for a theory to explain spiral structure in galaxies harmoniz- 
ing permanent structures in the disk with its differential rotation — thus avoiding the "winding dilemma" — and the 
Send offprint requests to: Fernando de Pablos 
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formulation of Lin (1968) in what would become the spiral density-wave theory, our conceptions about the origin of 
spiral structure and its relation with star formation have evolved, although not so dramatically as might be expected. 
Some points seem reasonably established, such as the presence of density waves in the disk (Elmegreen & Elmegreen 
1984), while others remain controversial. Can the arms propagate through corotation (Patsis et al.l994, and references 
therein)? Do density waves enhance star formation in the arms? Regarding this point, Elmegreen & Elmegreen (1984) 
found that grand-design spirals (those with continuous, well defined arms) have stronger density waves than flocculents 
(those with patchy and less well defined arms) , but that the arms of grand-design galaxies were nearly purely density 
wave enhancements. Moreover, from the analysis of several star formation tracers, Elmegreen & Elmegreen (1986) 
found that the star formation rate per unit area is virtually the same in flocculent that in grand-design galaxies, con- 
cluding that density waves do not significantly enhance star formation in the arms with respect to that in the interarm 
disk. Prom this and similar studies, such as that of McCall & Schmidt (1986), it seemed that the mechanism accounting 
for star formation in the arms would probably be self-propagating star formation as that proposed by Gerola & Seiden 
(1978), instead of a shock triggered when the material of the disk in its differential rotation encounters the density 
wave at supersonic speeds. This line of thought was reinforced by the findings of Lubow, Balbus & Cowie (1986), who 
showed that when taking into account the viscosity of the gas and/or gas gravity the gas does not shock. However, 
Cepa & Beckman (1990a), studied the efficiency of the massive star formation (say the massive star formation rate per 
unit mass of gas) in the arms with respect to that in the interarm disk (the so-called relative massive star formation 
efEciency) as a function of the galactocentric radius for some grand-design spirals and an intermediate-arm type one 
(M 33), finding that the spiral arms of the grand-design galaxies studied were more efficient at forming stars than the 
interarm disk by a factor ranging from 2 up to 15 or more, showing a pattern of peaks and troughs coinciding with the 
position of the resonances and corotation radius predicted by the density- wave theory, when a pattern speed is fitted to 
the curves Op = 0±K;/m, where Op is the pattern speed, is the rotation curve, k is the epicyclic frequency and m is 
the number of arms. This, and similar results that followed (Cepa & Beckman 1990c), led to the conclusion that some 
mechanism is bound to be forming stars in the spiral arms in a diff'erent way than in the interarm disk and consistently 
with the spiral density-wave predictions. Moreover, this behavior was not observed in the case of the intermediate-arm 
type galaxy studied (M 33), indicating that perhaps in at least some intermediate-arm type spirals the mechanism 
forming stars in the arms could be a self-propagating one, operating in the arms and in the interarm disk. But, even 
in this case, some spiral density wave has to be present in the disk to allow a self-stocastic star formation mechanism 
to form stable spiral arms, since three-dimensional simulations of this mechanism do not form stable spiral arms by 
themselves (Statler, Comins & Smith 1983). 

Some of these results are consistent with those observed by other authors using diflferent techniques (Tacconi & 
Young 1990 in NGC 6946 ; Wilson & Scoville 1991 in M 33). Although this could seem to be in contradiction with the 
studies of Elmegreen & Elmegreen (1986) and Lubow, et al. (1986), it was pointed out by Cepa & Beckman (1990b), 
on the basis of Elmegreen & Elmegreen's (1986) results, that if the star formation rate is not averaged over the entire 
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disk and only those parts of the disk where the star formation is taking place are taken into account, the star formation 
rate per unit area of star forming zones could be as much as 6 times larger in grand-design than in flocculent spirals. 
Also, star formation could be enhanced in the arms by mechanisms other than a large-scale shock (Elmegreen 1992, 
1994), and hence the expression "star formation front" would, in general, be more accurate than "shock front". 

However, the studies of relative arm/interarm star formation efficiencies mentioned above require spatially resolved 
measures of the Hi and CO over the whole disk of a galaxy, which limits the rate of data adquisition, thereby hampering 
statistically significant results unless large amounts of radio observing time are employed. There are also additional 
observational difficulties: the resolution required for this kind of work is critical and should be of the order of a star 
forming region, i.e., ~0.5 kpc (Cepa & Beckman 1990a), otherwise the arm/interarm CO contrast is underestimated. 
This resolution is difficult to achieve and slow to map in CO unless interferometric techniques are employed, in which 
case, and given the size of nearby galaxies, a mosaic of observations is needed (again slow to map). It could also be 
difficult to detect interarm CO using interferometers. Other critical points are that the conversion factor from CO 
to H2 (Garcia-Burillo & Guelin 1990), the initial mass function (see, for example, Shu, Adams & Lizano 1987, and 
references therein), and the physical properties of Hll regions, could be different in the arms from those in the interarm 
disk, and that the optical measures (to get the star formation rate) could be significantly affected by extinction. All 
these effects may cause severe alterations, either increasing or lowering the measures of arm/interarm efficiency ratios 
by a factor difficult to estimate, depending on the galaxy and the galaetoeentric radius within a given galaxy. 

The methods proposed in the present work allow the measurement of relative arm/interarm (or in general from 
one zone of the disk with respect to another) star formation rates and star formation efficiencies as a function of the 
galaetoeentric radius using optical broad- and narrow-band (Ha and H/3) imaging observations only. It is possible then 
not only to study these parameters for a larger sample of galaxies in a more efficient way, since no radio measurements 
are needed, but also, from the narrow-band data, to correct for extinction and to diagnose possible differences between 
the initial mass function in the arms and the interarm zone, thus avoiding most of the uncertainties mentioned in the 
previous paragraph. 

The method is described in § 2. Section 3 describes the observations, data reduction, flux calibration and the 
procedure followed to apply the method to the observed spirals, while the results are discussed in § 4. Finally, a 

summary is given in § 5. In this paper we present the method and an example of the application to two grand-design 
spirals. In Paper II we will present the results obtained for a larger sample of objects of different arm types. 



2. Method 
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2.1. Relative star formation rates 

We consider the total luminosity of a massive star forming zone (i.e., an Hii region) in the disk of a galaxy as consisting 
of the luminosity due to the massive star formation plus the luminosity of the underlying older disk population. The 
total luminosity at a given A of a massive star forming zone in the interarm disk can then be written as 

^lARM ^ ^DISK ^ ^DSF 

where L^^^^ is the luminosity, at a given A, of the underlying older disk, and L^^^ , the luminosity due to the massive 
star formation in the interarm disk. Also, the total luminosity of a massive star forming zone in the arm is, 

where L^^^ is the luminosity due to the massive star formation, and is the stellar density contrast (in the arm 
with respect to the interarm disk). All these quantities depend on the galactocentric radius although the dependence 
is not explicit to simplify the equations. 

The luminosity due to recent star formation can be expressed as a function of the star formation rate (SFR) 
the initial mass function (IMF), and the luminosity, at a given A, of the stars of masses between m and m + dm, 



Lx^Tp (j){m)Lxim)dm^ij<^x (3) 

J mi 

Then represent the luminosity per unit mass emitted at A of the stars formed, and (|l|) and (|^) can be written, 

LlARM ^ ^DISK ^ ^^DSF <j,DSF (4) 

and 

L^I^^ = A. (L^ISK ^ ^ASF ^ASF) (5) 

ip-^^^ and ^p^^^ are the SFRs in the arm and in the interarm disk, ^x depends on the IMF and the luminosity of the 
stars in every range of masses at the given passband A, and we will consider that the only possible differences between 
$A in the arm and the interarm region will be caused by differences between the IMFs in the arm and the interarm 
region due to a biased initial mass function. 

Using these equations, the relative arm over interarm star formation rate (RSFR) may be expressed, 

A ,/>ASF rARM A r DISK rhDSF 
^ l/)DSF 2^IARM _ T^DISK ^ASF 

L^™, i^i^ARM j^j^j ^DiSK g^j^g observable quantities that can be obtained performing aperture photometry on a spiral 
arm, a massive star forming zone in the interarm disk, and a zone of the interarm disk with little or no star formation, 
respectively. The passband A has to be sensitive to star formation, i.e., it has to allow the clear distinction between star 
forming zones and non-star forming zones. In this work we will use the Johnson B band, which fulfils this requirement 
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and is, at the same time, easy to observe (for example, the Johnson U band requires much longer integration times 
to achieve the same signal-to-noise ratio). Also, is reasonably free of emission lines: only H/3 lies on the red edge of 
the filter spectral response curve (Allen 1976) at ~50% of the peak response, or even less, given a redshifted systemic 
velocity. 

The stellar density contrast can be estimated from the ratio of arm to interarm disk luminosities of non-star forming 
zones in a passband less sensitive to star formation, i.e., that approximately traces stellar density variations. We have 
used the / band, because it is less affected by extinction and is reasonably free of strong emission lines: [O II] at AA 
7325 A lies on the blue edge of the filter at ^65 % of the peak response (Benn & Cooper 1987), and [S III] at AA 9069 
A falls on the red edge at ~25 % of the peak response (Benn & Cooper 1987), or even less given a redshifted systemic 
velocity. Then, 

jARM 

where Kj is a factor to correct for the contribution to / band from newly formed stars in the arm. From Schweizer 
(1976) it results that Kj ^ 0.85, in agreement with the observational result of Elmegreen & Elmegreen (1984). 
Finally, using (0) and applying (^) to Johnson B band, (^) can be written, 

^ARM jDISK _ ^^jARM^DISK ^ 
RSFR — jDISK^^IARM _ ^DISK^ ~ 

In general, xb = ^^^^/^^^^ > 1. If the IMF is the same over all the disk, xb = 1- If the star formation is biased 
towards a larger fraction of massive stars in the arm with respect to the interarm disk, xb > 1- Moreover, in this case 
the value of xa depends on the passband used: the effect of biased star formation (BSF) would be more noticeable at 
shorter wavelengths. 

All the quantities on the right hand side of (||) are luminosities which can be measured directly, except xb, which 
is not an observable in the present work. Then, 

and RSFR™'''*^ is the RSFR if no BSF is present, otherwise the factor xb has to be taken into account. Also, RSFR™''^'' 
gives an estimate of the relative star formation rate for stars of intermediate masses, i.e., the ones that contribute 
more to the B-band luminosity. 

2.2. Relative star formation efficiency 

We define the relative (arm with respect to the interarm disk) star formation efficiency (RSFE) caused by density 
waves as 

RSFE^^ (10) 
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which represents the efBciency of density waves inducing star formation in the arms. From (^): RSFE= ^jArm^^disk_ 
As in the case of RSFR, 

Henceforward, the RSFR™°'^'' and RSFE™''^'' wiU be termed intermediate-mass star formation rate and efficiency, 
respectively. 

2.3. Relative massive star formation rates 

The Ha luminosity (Lhq) of a star forming zone can be related with the number of ionizing photons (A^uv), assuming 
case B conditions and spherical Hii regions via (Osterbrock 1989): 

Lna = /iz^HaiVuve— (12) 

where ly^a is the frequency of the Balmer a line, e measures the efficiency of the gas absorbing photons, and and 
as depend on the temperature (and hence on the position of the Hii region), and their values are tabulated for case 
B recombination (Osterbrock 1989). 

By assuming a certain IMF it is possible to relate the Ha luminosity with the star formation rate: 



a 



.off rm 
Ha 



Laa — hvuai'^ — — / 4)(rn)Nuy{m) dm (13) 

J ra-B 

where me is the lower mass limit for a star to emit photons with a frequency above the Lyman limit, and Nuy^m) is 
the number of ionizing photons emitted by a star of mass to. This relation will, in general, depend on the position in 
the disk (due, for example, to different metallicities). Also, 

N^y{m)= r^^^d,. (14) 

where L(to, ly) is the luminosity of a star of mass to at frequency v, and Vo is the frequency of the Lyman limit. 
Then, if 

,1 ARM I ARM ,efT ARM 

n = ^ d 111 

^ ,ARM^ARM,off lARM ^ ' 

and 

JruB ^ ^uv dm 

the relative arm/interarm star formation rate inferred from massive stars can be written, using the same notation as 
in the previous subsection, 

^^.0ASF Ha^^^ 1 

RSFRoB = = Q n^iARM — (17) 
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where Ha^^*^ and Ha^^^*^ are, respectively, the Ha luminosities of a star forming region in an arm and in the interarm 
disk, respectively, at a similar galactocentric distance R (because quantities in (|l7| ) and ( |l5| ) depend on the position 
in the disk). 

It is difficult to estimate the ratio Q of (p^). Let us consider first the ratio between interarm and arm zone of the 
term a^^/oB- This ratio depends on the temperature. The densest Hii regions could be situated in the arms, and 
densest Hii regions are likely to be hotter (McCall, Rybski & Shields 1985). If we assume that in the most extreme case 
Hii regions would have a temperature of 15000 K in the arms and 5000 K in the interarm regions, from Oster brock 
(1989) it turns out that this ratio takes a value of 1.1. However, the metallicity of Hii regions in the arms might be 
higher, and this effect would tend to lower the temperature (Shields 1990, and references therein). Also, the mean 
temperature of Hii regions in the arms may well be similar to that in the interarm zone, and then the ratio would be 
of order unity. We can then consider that 

IARM,cff ARM 
"Ha "b 

Concerning the ratio between the gas efficiencies in absorbing photons in arm and interarm zones, McCall et al. 
(1985) pointed out that the Hii regions of their sample were ionization bounded, so that in this case e would be 1.0. 
Moreover, there is prima facie evidence that the Hii regions in the arms and in the interarm zone do not differ in 
their boundary conditions (Cepa & Beckman 1989; Cepa & Beckman 1990d), thus leading to the conclusion that it is 
likely that giARM/gARM ^ ^ q However, the regions of the sample of McCall et al. were mainly located in the arms, 
and it may be argued that, at least in some cases, Hii regions in the interarm zone are density bounded while Hii 
regions in the arms are ionization bounded. In this case the ratio jIARm^^arm -^^ould be less than unity. We can do a 
first-order estimate of this ratio, taking e as proportional to the number of atoms divided by the number of ionizing 
photons in the Hii region. The column density of atoms that can be ionized in a cloud is proportional to its diameter 
and density, and the number of ionizing photons is proportional to the SFR and IMF. The SFR is itself proportional 
to the product of the gas density and the star formation efficiency (SFE). If the SFE is larger in the arms and/or 
the IMF is different in arms and interarm regions, in the sense of favoring the formation of more massive stars in the 
arms, then giARM^^ARM -> ^iarm^^arm^ where (p represents the diameter of the Hii region. For example, from Cepa 
& Beckman (1990d), this ratio is 0.75 for NGC 4321, and from Knapen et al. (1993) this ratio is 0.72 for NGC 6814. 
Then, in general, 

lARM 

0.7 < < 1.0 (19) 

and finally, 



0.7 < g < 1.1 



(20) 
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This factor has to be taken into account when evaluating relative arm/interarm star formation efficiencies using 
Ha luminosities. In the present work, 

RSFRJ^j;^"^ = RSFRoB ^ (21) 
and 

RSFEgl^^" = (22) 

Note that, except for the factor Q, RSFRob'" and RSFEqb''' are an estimate of the relative (arm over interarm disk) 
star formation rate and efficiency, respectively, of massive stars, and so will be termed henceforward when referring to 
the measured quantities. 

2.4- Biased star formation 

The ratio between RSFR and RSFRqb has to be unity. However, this does not apply to the ratio of the measured 

RSFRs, 

RSFR^ 
rtor riQg 

where x = Xb/xho- 

If there is no biased star formation, x — ^- Otherwise, x would be less than 1.0, because Ha samples higher stellar 



masses than B. Then the ratio (23) could provide an observational test for the presence of biased star formation. From 
(©, 

0.7 V < < 1.1 Y (24) 

tXDf Jaq3 

so it is possible to affirm that there is evidence for biased star formation in the arms with respect to the interarm disk 
if RSFR^^'^VRSFRJ^^'^" < 0.7. 

3. Observations, data reduction and flux calibration 

3.1. Observations 

In this paper we present some results derived from the application of the method previously presented to two grand- 
design spiral galaxies, NGC 4321 and NGC 4254 (Table 1). 

Images were taken for both galaxies in broad (B and /) and narrow bands (Ha and H/3) . For narrow-band images 
we used interference filters, including the corresponding continuum filters to the lines Ha and H/3. To derive the filter 
transmission in the Ha and H/3 lines, both the redshift of the galaxy as well as the shifts in wavelength in the filter 
due to temperature changes (Peletier 1994) were taken into account and, when applicable (as in the case of the prime 
focus of the Isaac Newton Telescope), the shift in central wavelength due to the non-converging beam (Peletier 1994, 
see Tables 2 and 3). 
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The images of NGC 4321 in the B, I, Ha and Ha continuum filters were taken in 1994 March with a 1024x1024 
pixel Thompson CCD camera attached to the cassegrain focus of the 0.8-m IAC-80 telescope at the Teide Observatory 
(Tenerife, Spain), the scale obtained with this configuration being 0.43 arcsec pix~^. The H/3 and H/J continuum images 
were obtained in 1995 June at the Cassegrain focus of the 1.5-m Danish Telescope at ESO (La Silla, Chile), using a 
1024x1024 pixel Tektronics CCD, giving 0.38 arcsec pix~^. The images of NGC 4254 were taken in 1995 April with 
a 1024x1024 pixel Tektronics CCD camera attached to the prime focus of the 2.5-m Isaac Newton Telescope (INT) 
at the Roquc dc los Muchachos Observatory (La Palma, Spain). In this case the scale obtained was 0.59 arcsec pix~^. 
Landolt catalogue stars (Landolt 1973) were used for the broad-band calibration. For the narrow band wo used the 
spectrophotometric standard stars of Okc & Gunn (1983) and Stone (1977). Atmospheric extinction was evaluated in 
the broad-band filters via the observation of Landolt (1973) standards at different airmasses each night or for the INT 
data using the V extinction given by the Carlsberg Automatic Meridian Circle. These extinctions allow the evaluation 
of extinctions in other filters assuming a grey atmosphere and a theoretical extinction curve versus wavelength (King 
1985). In Table 4 we present the total integration times, and in Table 5 the limiting fluxes and magnitudes attained 
in the observations. To facilitate the continuum substraction and to correct for underlying absorption lines (described 
in § 3.3 and § 3.4), both the on-line and the continuum images were flux calibrated. 

3.2. Data reduction 

Data reduction was performed using the IRAF package. The bias was subtracted using the overscanned part of the 
chip. Several flatfields were taken for each filter each night and were later averaged to obtain a master fiatfield for 
each filter and each night. The sky background was determined for each image by averaging and substracting regions 
free of objects. Afterwards, in the case of the observations taken with the IAC-80 or the ESO Danish Telescope, where 
several 1800 s exposures were taken for each filter (except for the continuum H/3 at the ESO Danish Telescope where 
the exposure time was 1200 s per image), the off-sets among the different images were determined by fitting Gaussians 
to several field stars in each frame, the images were aligned to a fraction of a pixel and then averaged using a sigma 
clipping algorithm to get rid of cosmic ray events. The single exposures obtained with the INT were interactively 
cleaned of cosmic rays. 

3.3. Flux calibration and continuum subtraction 

Both on-line and continuum images were flux calibrated using a procedure similar to that described in Earth et al. 
(1994). The flux calibrated continuum image was then shifted to a fraction of a pixel to match the on-line image, 
scaled by different factors around 1.00, in steps of 0.05, and substracted from the on-line image. From the different 
resulting images, the one that showed values closer to zero in regions of the interarm disk free from Hll regions, and 
no conspicuous negative values in the bulge zone was selected. The scaling factors found in this way were 1.00 for 
Ha images and 0.90 for H/3 images. This procedure for the continuum substraction might eliminate diffuse emission 
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coming from the interarm disk, but this contribution is irrelevant for the present work (see di Serego AHghieri 1987 
for an overview of on-hne imaging and its pitfalls). Figures 1 and 2 show the resulting Ha images. 

3.4- Application of the method 

To apply the method outlined here to the data, it is necessary to evaluate luminosities through three different sets of 
apertures for each galaxy: 

1. The apertures to measure arm luminosities have been taken as contiguous, except when there are abundant arm 
Hll regions, where the apertures overlap slightly, and when very few or no arm Hil regions are apparent, where the 
apertures undersample the arm. These arm apertures are located in the most intense part of the Hq image of each 
arm. When the Ha image presented a gap in the arms, the S-band image was used as a guide to follow the arm 
in the Ha image. 

2. The apertures to measure luminosities of star forming (Hii) regions on the disk have been defined using the Ha 
image. Zones of doubtful assignation (to an arm or to the disk), too near to dust lanes or to the tips of a bar, have 
been avoided. 

3. Finally, aperture positions for zones of the disk without star formation have been selected avoiding: (i) disk star 
forming zones (the Ha image has been used for this selection), (ii) emission coming from the spiral arms (using the 
B and / images) , and (iii) dust lanes (using the B and / images) . The latter is necessary since it is not possible to 
correct these points for internal extinction, only for Galactic extinction. Except for the previous stated exceptions, 
zones equidistant from the spiral arms and avoiding bars have been used. 

The size of the apertures employed is fixed for each galaxy and is large enough to get isolated Ha emission patches 
of the disk fully inside them. This implies diameters of 8" for NGC 4321 and NGC 4254 (~0.8 kpc assuming a distance 
of 20 Mpc). 

Afterwards, aperture photometry for the whole set of apertures and all the images (except Ha and H/3 and their 
corresponding continua in the case of disk regions without star formation) was performed. The coordinates of the 
apertures with respect to the center of the galaxy, defined as the peak of maximum emission in the / band, were then 
deprojected using the position angles and inclinations given in Table 1. 

The luminosities thus obtained were corrected for extinction and underlying absorption (in the case of H/3) following 
the procedures described below. 

3.4.1. Extinction correction 

Since extinction in the disk can be quite large, especially in the spiral arms, the relative star formation rates must 
necessarily be corrected for internal extinction to be able to apply the method described in § 2 with a certain confidence. 
For instance, extinction can either increase or decrease the RSFR, the RSFE and the RSFEqb, depending on the zones 
considered in the arms and in the interarm disk, while usually RSFRqbj and A^, will be lowered by extinction (assuming 
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that the extinction is larger in the arms than in the interarm disk). In the present work, we have used the ratio of 
the Balmer emission lines Ha and H/3 to do a first-order correction for galactic and extragalactic extinction, using the 
reddening law of Whittet (1992), assuming case B with an electron density of 100 cm~^ and a temperature of 10"* K, 
from which we obtain 



which allows us to derive extinctions for our photometry. 

This reddening law is not substantially different from that of Schild (1977). Moreover, recent results (see Davies 
& Burstein 1995) show that it can be applied reliably to external galaxies, at least in the visible and NIR bands. 

However, this method can be applied only to Hii regions. So, with this method, the luminosities with a superscript 
DISK (§ 2) cannot be corrected for internal extinction, only for extinction due to our Galaxy and internal extinction 
for correction to face-on, using Ab from de Vaucouleurs, de Vaucouleurs & Corwin (1976). In these circumstances, 
RSFEgg'" is lower, the measured is larger and if RSFE"^'=''" < 1 (see Appendix) RSFR'"'^'''^ and RSFE'"'=^" are 
larger than the same fully extinction-corrected quantities. Also, if RSFE'"^''^ < 1, the ratio given by equation (|2^) is 
larger than if it were fully extinction corrected, and then it is still true that if ( p3|) is lower than 0.7, there is evidence 
for biased star formation. 

3.4.2. Line-absorption correction 

Apart from extinction, there are several effects that may influence the observed fluxes, namely the absorption lines in 
the calibration stars and the underlying stellar absorption of the studied region. The first effect is negligible because 
the equivalent widths of the absorption lines are ~2 A (McCall et al. 1985) i.e., small in comparison to the filter 
widths. Hence, the absorbed flux is always lower than 4% of the continuum flux for a filter of 50 A FWHM or more, 
such as those used in this work. The second effect may be more significant reaching reductions in the H/3 Balmer 
line equivalent widths of 50% (McCall et al. 1985). This effect leads to an overestimate of the absorption. For this 
correction we used the method given by McCall et al. (1985), which allows us to correct for the underlying absorption 
effect using the equivalent widths of the Ha and H/3 lines. First, a mask was derived with pixel values equal to unity 
when the pixels of the same coordinates in the Ha line image had values greater than or equal to three times the 
r.m.s. of the background, and pixel values otherwise equal to zero. This mask was multiplied by the continuum images 
in order to obtain continuum the contribution of emission-line regions only, not of regions which are pure continuum. 
This allows to evaluate equivalent widths. Since the internal (non-atmospheric) extinction is different in the line and in 
the continuum due to the different central wavelengths of the filters used, it is necessary to apply an iterative process 
to correct for underlying absorption and to derive extinctions: initially the Ha and H/3 equivalent widths uncorrected 
for extinction are calculated, from those we evaluate the absorption-corrected H/3 ffuxes, before deriving the internal 
extinction. This internal extinction is used to correct the Ha and H/3 line and continuum fiuxes. Using these first-order 



Av = 2.6 In 



(-FHa/-pHff)Qbs 



(25) 



2.86 
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extinction-corrected fluxes, we calculate the equivalent widths again, and the H/3 flux is later corrected for underlying 
absorption, which allows the determination of the extinction, etc. This procedure is repeated until a convergent value 
is obtained, which takes only few iterations. 

Another effect that might lead to errors in line intensities is contamination produced by the [N II] A 6584 A line. 
However this effect will not affect the arm/interarm ratios, except in the case that the variation of the metallicity 
between the arm and interarm is significant, since arm/interarm ratios are evaluated at similar galactocentric distances, 
thus avoiding the effect of possible metallicity gradients in the disk. 

Finally, the deprojected, extinetion-corrected luminosities, for each arm, and interarm disk, with and without star 
formation, at similar deprojected galactocentric distances, were used to apply the formulae of § 2, to determine arm 
amplitudes, intermediate-mass star formation rates, massive star formation rates, the corresponding star formation 
efficiencies, and look for possible biased star formation in spiral arms with respect to the interarm disk as a function 
of radius. 

4. Results 

4.1. NGC 4321 

NGC 4321 is one of the most conspicuous spiral galaxies in the Virgo cluster. It has two impressive, very well defined 
spiral arms (of arm class 12 according to Elmegreen & Elmegrecn 1987). The northern arm extends in the east-north- 
west direction and the southern arm in the west-south-east direction. 

4.1.1. Triggered star formation 

Figures 3 and 4 represent the relative (arm/interarm) intermediate-mass star formation rate, and the relative massive 
star formation rate versus radius for the southern and northern arm, respectively. In both arms the qualitative behavior 
of both quantities as a function of radius is noticeably similar for each arm. In the southern arm there is a descent 
below unity of both RSFRs from 105" to 140". In the northern arm there is a dip at 110" and a descent below unity 
at 140". In these zones, the stellar arm density contrasts are also lower (Fig. 5). In the zones where the RSFRs are 
less than unity, the interarm disk is forming more stars per unit time than the arm, and can be identified with parts 
of the arms with almost no Hii regions, as can be seen in the Ha image (Fig. 1). These dips and descents can be due 
to the presence of corotation or to resonances. In the modal theory due to Lin and collaborators, oscillations in the 
amplitude of the wave are related to the oscillations of the radial eigenfunctions (see for example Lin & Lau 1979). In 
the stellar theory for non-linear spirals (Contopoulos & Grosb0l 1986, 1988) there is evidence for stellar orbit crowding 
between resonances, and it is possible to find some minima in the arm amplitudes placed at different resonances. We 
assume that corotation is located at ~110-120", where dips in the RSFRs, and stellar density contrasts can be seen 
in both arms. In this case, the dip at 140" corresponds to the 4:1 Outer Lindblad Resonance (04:1), while the dip 
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at 60" coincides with the 4:1 Inner Lindblad Resonance, using the rotation curve published by Sempere et al. (1995). 
The corotation radius determined in this way is consistent with that of Sempere et al. (1995), who found that it is 
located from 85" to 115", using a test based on the Hi residual velocity field, and that of Elmegreen, Elmegreen & 
Seiden (1989), at ~120", using morphological arguments. 

Figure 6 shows the relative intermediate-mass star formation efBciency versus radius for both arms of NGC 4321. 
The arms are less efficient, or as efficient at forming intermediate-mass stars than the interarm (then the reasoning 
of Appendix A is valid, implying that the RSFE of massive stars is lower, and A^, RSFR and RSFE of intermediate- 
mass stars are larger than the same quantities using extinction-corrected disk non-star forming zone luminosities), 
except beyond 150 arcsec, where measurement uncertainties are quite large. In contrast. Fig. 7 shows the relative 
massive star formation efficiency versus radius for both arms of NGC 4321. In some zones, the southern arm is more 
efficient at forming massive stars than the interarm disk by a factor up to -^2.5-3.5, indicating that density waves 
are triggering massive star formation, a result which was already pointed out by Cepa & Beckman (1990c) using gas 
density contrasts. However, in the northern arm, the RSFEqb''* is, in general, less than one. Del Rfo (1995) and del 
Rfo & Cepa (1997), using azimuthal photometric profiles in Johnson B and / bands, found that the density wave is 
quite inefficient at forming stars in the northern arm, except in a localized region outside corotation, which can be 
identified with the region of high RSFEgg''' at ~150" in Fig. 7. Nevertheless, arm density contrasts are in general 
larger in the northern arm than in the southern arm (Fig. 5), the reason probably being due to the physical conditions 
of the molecular gas. Also, if there is no triggering of massive star formation, a linear relation between the RSFRgg***, 
and the arm density contrasts should be expected, and vice versa: triggering should show up as a non-linear relation 
between arm density contrast amplitude and relative massive star formation. From Fig. 8 it appears that there is no 
evidence for correlation between the RSFRob*^ and the southern arm amplitude, again leading to the presence of some 
sort of triggering mechanism. However, there is a remarkable linear relation between the points where no triggering of 
massive stars is present (with a RSFEqb*'* less than one), which are situated around the corotation radius (Pearson 
correlation coefficient rp = 0.993 for n = 5, giving a confidence level of 99.9 % for this correlation). In the northern 
arm, avoiding the points being beyond 140" (where there is some triggering of the star formation), we observe a 
linear relation between RSFRq^*** and the north arm density constrast (Fig. 9), with a Pearson correlation coefficient 
rp = 0.891 for n = 9, at a 99.8 % confidence level for this particular correlation. This linear correlation points to the 
absence of massive star formation triggering in the northern arm, since larger arm density contrast corresponds to 
proportionally larger RSFRqb^''. However, there is evidence for triggering of the massive star formation beyond 140". 
This difference in behavior between both arms might be due to the interaction of NGC 4321 with NGC 4322. 

4.1.2. Biased star formation 

Prom the previous section, it turns out that the density wave system in NGC 4321 is triggering the formation of 
massive stars, while intermediate-mass stars are formed in larger fractions in the interarm disk. This result suggests 
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that there might be biased star formation in the spiral arms with respect to the rest of the disk. As pointed out in 
§ 2.4, values of RSFR™®'*^/RSFRg|*'* lower than 0.7 constitute strong evidence for a different IMF in the arms with 
respect to the interarm disk, or, in other words, biased star formation (BSF). In Fig. 10 we have plotted this ratio as a 
function of the galactocentric radius for the northern and the southern arms. In the northern arm, the ratio of relative 
measured star formation rates is only systematically lower than 0.7 from 130" to 150", and in 60" and 90". It would 
then seem that there is BSF, and that it is mainly situated in the zone where there is triggered massive star formation. 
In the southern arm the ratio is less than 0.7 except in a zone from 100" to 130". Again, there is BSF, mainly in the 
zone where massive star formation is triggered, and not near corotation (where no triggering is expected). This is a 
strong suggestion that density waves cause a change in the IMF in the arms. This is in some agreement with Shu, 
Lizano & Adams (1986) who proposed that regions of higher (lower) star formation efficiency corresponds to regions of 
high (low)-mass star formation. A possible reason might be (Zinnecker 1987) that the giant molecular clouds — located 
in the spiral arms — form massive stars earlier than the standard molecular clouds, the less massive stars not having a 
chance to form. Also this result is in agreement with that of Scoville, Sanders & Clemens (1986) who proposed that 
OB stars form as a result of cloud-cloud collisions, which happen more frequently in the arms. There is some other 
evidence for the presence of BSF in spiral galaxies: Miller & Scalo (1978) conclude that in OB associations in our own 
Galaxy, which are tracers of the spiral arms, only stars with M > 2 — 5 Mq are formed. Beyond our Galaxy Rieke et 
al. (1980) found a burst in the center of M82 where they estimate a value of nii ~ 3.5 M© for the lower mass limit 
(see the review by Shu et al. 1987 for more examples). 

4.2. NGC 4254 

lye et al. (1982) found a strong m =1 component in the Fourier transforms, superimposed on the m = 2 component, 
giving the observed m = 3 spiral structure in NGC 4254. Recent studies made by Phookun, Vogel & Mundy (1993) 
ascribe the cause of the strong m = 1 component to the presence of an extragalactic cloud which merges with the 

galaxy. 

The optical images clearly show three arms, one of them appearing weaker in Ha than the others. We have studied 
the star formation in the strongest arms: the northern arm, which extends in the west-north-east direction, and the 
southern arm, which extends in the east-south- west direction (Fig. 2). 

4.2.1. Triggered star formation 

In Figs. 11 and 12 we have plotted RSFR™'"*'* and RSFRqb'^ versus radius for the south and the north arms, respec- 
tively. In both arms there is a common minimum of both quantities between 85" and 95", which also appears in the 
arm amplitudes and in RSFEgg^'', as can be seen in Figs. 13 and 14, respectively. From this behavior we should expect 
that the corotation will be between 85" and 95", a result consistent with findings reported by Elmegreen, Elmegreen 
& Montenegro (1992), who placed the corotation at ^ 87". Using this corotation radius, from Fig. 14 it turns out 
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that RSFEqb'^ in the northern arm is greater than unity before and after corotation, and less than unity near or at 
corotation, as in the southern arm of NGC 4321. This arm of NGC 4254 contains a region of strong star formation 
at 115" (Fig. 2), which has not been included since it is quite possible that it was not triggered by the density wave. 
In fact, the massive star formation efficiency in that zone is an order of magnitude greater than in any other zone of 
star formation in this galaxy, and it coincides approximately with the area where the extragalactic cloud observed by 
Phookun et al. (1993) merges with the galaxy, this merger being perhaps the cause of the strong star formation in the 
zone. In the southern arm (Fig. 14) there is also some triggering of massive star formation mainly after corotation. 

In general, there is no evident correlation between RSFRgg'^^ and the stellar density contrast in any arm (Figs. 
15 and 16). This fact supports the hypothesis that there is triggering of massive star formation in the spiral arms 
with respect to the interarm disk. However, there is a good linear relation for the points in the southern arm having a 
RSFEol"*^ clearly less than one (Fig. 15), i.e., zones where there is no triggering of massive stars (Pearson correlation 
coefficient rp = 0.994 for n = 5, giving a confidence level of 99.9 % for this correlation). 

As in the case of NGC 4321, the measured RSFEs (Fig. 17) are less than or of the order of unity except in some 
points of the southern arm at ~30, 40, and 60 arcsec. The density wave system, then, is triggering massive star 
formation in the arms, but intermediate-mass stars are formed preferentially outside the arms, suggesting again the 
possible presence of biased star formation in the arms. 

4.2.2. Biased star formation 

Figure 18 represents the ratio RSFR^^^^/RSFRgl^^ as a function of the galactocentric radius for the south and 
the north arm. This value is lower than 0.7 in both arms beyond corotation, implying strong evidence for BSF after 
corotation, although a BSF before corotation for the southern arm (depending on the values of Q) cannot be ruled out. 
The lowest value corresponds to 120" in the southern arm, correlated with the peak in RSFEgg*^ at this distance (Fig. 
14). Although there seems to be no correspondence between the peak in RSFEgg*^ at --^SO" observed in the northern 
arm and a dip in RSFR™''*'*/RSFRqb^'^, the observed change of behavior in this ratio before and after corotation, 
mainly in the southern arm, suggests that the density wave is responsible for the observed BSF. 

5. Summary 

We have developed a procedure for deriving relative (arm/inter arm) star formation rates, arm density contrasts, and 
relative star formation efficiencies as a function of the galactocentric radius in spiral galaxies, using broad-band (blue 
and near infrared) and continuum-subtracted Ha and H/3 imaging. We have also developed a procedure for assessing 
whether there is biased star formation in the arms with respect the rest of the disk, i.e., whether the spiral arms have 
a different initial mass function from that of the interarm disk, as a function of radius. 

We have applied the previous procedures to the spiral galaxies NGC 4321 and NGC 4254, finding that density 
waves greatly influence the star formation and initial mass function in the spiral arms. 
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Prom the relative star formation rates and arm amplitudes as a function of the radius, it results that the corotation 
radius is located at 110-120" in NGC 4321 and at 85-95" in NGC 4254. 

Prom the relative intermediate- mass star formation efficiencies and massive star formation efficiencies as a function 
of radius, it is found that density waves trigger massive star formation in both galaxies, although not with the same 
strength in the different arms of the same galaxy. However, they do not trigger intermediate-mass star formation. 
Also, triggering is not present near the corotation radius where, in the case of NGC 4321, a linear relation between 
relative massive star formation rate and arm amplitude can be observed. This linear relation is also apparent in the 
zones outside corotation where the relative massive star formation efficiency is less than unity, i.e., where no triggering 
of massive star formation by the density wave is taking place. 

Near corotation (and, in the case of NGC 4321, near the Inner and Outer Lindblad 4:1 Resonances), the interarm 
disk presents larger star formation rates and star formation efiiciencies than the spiral arms, as if the star formation 
were inhibited in the spiral arms near corotation. In fact, a break in the arms or a clearly lower number of Hll regions 
in this zone can be seen in the Ha images. 

Intermediate-mass stars form more efficiently in the interarm disk, and the IMF is different in the arms from that 
in the disk, favoring the formation of a larger fraction of massive stars in the arms. This biased star formation seems 
closely related and probably caused by the density-wave system, since it is not present where there is no triggering of 
massive star formation (relative massive star formation efficiencies less than unity), in the case of NGC 4321, or it is 
present after the corotation radius and not before, as in NGC 4254. 

The observed correlation between biased star formation and higher massive star formation efficiency is consistent 
with the theoretical work of Shu et al. (1986). 

Pinally, we can conclude that, for the grand-design galaxies studied, density waves trigger massive star formation, 
but do not trigger star formation of stars of all masses, on the contrary, the interarm disk is more efficient at forming 
intermediate-mass stars than the arms, and this is due to different IMFs in the arms from those in the rest of the disk. 

To be able to draw more general conclusions, it is necessary to apply the method described in this paper to a 
larger sample of objects, including intermediate-arm type spirals, where no triggering and consequently no biased star 
formation should be expected. 
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A. Appendix 

Differentiating (||) considering as variables /^isk ^disk^ ^^^^ dividing the result by (||), we obtain 

ARSFR _ -AA.gDiSK ^^^^disk A^disk 

RSFR " BARM _ A^BDiSK 5 arm _ a^^disk + ^iarm _ ^disk 

where 

AA, _ A/DISK 

A, ^ /DISK 

AfiDiSK g^j-^j ^jDiSK j^j^g negative quantities that represent the correction to apply to the extinction corrected S^isk 
and /DiSK^ respectively, to obtain the measured corresponding quantities. Then, AA^, > 0. 

However, since extinction is larger at shorter wavelengths, and the colors of disk non star forming regions are red, 



(A2) 



|^jDiSK|/jDiSK < |AsDiSK|/^DiSK^ ^^le second term of (O) is larger than the first 



If RSFE < 1, then RSFR < A^,, and, neglecting the first term of (|AlD, it results that ARSFR > 



Using a similar procedure, it can be demonstrated that if RSFE< 1, then ARSFE > 0. 
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Figure captions 

Fig. 1. Ha continuum-subtracted image of NGC 4321. North is top and east left. Axes represent offsets in arcseconds 
from the nucleus 

Fig. 2. Hq continuum-subtracted image of NGC 4254. North is top and east left. Axes represent offsets in arcseconds 
from the nucleus. The star formation region marked with an "X" sign has not been included in the present study 

Fig. 3. Relative massive star formation rate and relative intermediate-mass star formation rate as a function of the 
galactocentric radius for the southern arm of NGC 4321 

Fig. 4. Relative massive star formation rate and relative intermediate-mass star formation rate as a function of the 

galactocentric radius for the northern arm of NGC 4321 

Fig. 5. Arm amplitudes of NGC 4321 as a function of the galactocentric radius 

Fig. 6. Relative intermediate-mass star formation efficiency for the spiral arms of NGC 4321 as a function of the 
galactocentric radius 

Fig. 7. Relative massive star formation efficiency for the spiral arms of NGC 4321 as a function of the galactocentric 
radius 

Fig. 8. Relative massive star formation rate as a function of the arm amplitude for the southern arm of NGC 4321 

Fig. 9. Relative massive star formation rate as a function of the arm amplitude for the northern arm of NGC 4321 

Fig. 10. Ratio of the relative (arm/interarm) intermediate-mass star formation rate over that of massive stars as a 
function of the galactocentric radius for the northern (upper diagram) and southern (lower diagram) arms of NGC 
4321. Points below the continuous line (which indicates the 0.7 value) indicate the presence of biased star formation 

Fig. 11. Relative massive star formation rate and relative intermediate-mass star formation rate as a function of the 
galactocentric radius for the southern arm of NGC 4254 
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Fig. 12. Relative massive star formation rate and relative intermediate-mass star formation rate as a function of the 
galactocentric radius for the northern arm of NGC 4254 

Fig. 13. Arm amplitudes of NGC 4254 as a function of the galactocentric radius 

Fig. 14. Relative massive star formation efBciency for the spiral arms of NGC 4254 as a function of the galactocentric 
radius 

Fig. 15. Relative massive star formation rate as a function of the arm amplitude for the southern arm of NGC 4254 

Fig. 16. Relative massive star formation rate as a function of the arm amplitude for the northern arm of NGC 4254 

Fig. 17. Relative intermediate-mass star formation efficiency for the spiral arms of NGC 4254 as a function of the 
galactocentric radius 

Fig. 18. Ratio of the relative (arm/interarm) intermediate mass star formation rate over that of massive stars as a 

function of the galactocentric radius for the southern (upper diagram) and northern (lower diagram) arms of NGC 
4254. Points below the continuous line (which indicates the 0.7 value) indicate the presence of biased star formation 
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Table 1. Observational parameters of the galaxies observed 



Galaxy 


Hubble type 


Arm class 


P.A. 


i 


Db25 








(deg) 


(deg) 


(arcsec) 


NGC4321 


SAB(s)bc(l) 


12(2) 


151(3) 


31(3) 


6.9(1) 


NGC4254 


SA(s)c(l) 


9(2) 


58(1) 


27(1) 


5.4(1) 



(1) de Vaucouleurs, de Vaucouleurs & Corwin (1976) 

(2) Elmegreen & Elmegreen (1984) 

(3) Cepa et al. (1992) 
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Table 2. Characteristics of the narrow band filters for NGC 4321 



Filter 


Ac 


FWHM 


Peak response 




(A) 


(A) 




Hq 


6570 


95 


0.63 


H/3 


4897 


50 


0.78 


HaC 


6260 


133 


0.53 


H/3C 


4926 


53 


0.67 
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Table 3. Characteristics of the narrow band filters for NGC 4254 



Filter 


Ac 


FWHM 


Peak response 




(A) 


(A) 




Hq 


6607 


53 


0.54 


H/3 


4863 


97 


0.68 


HaC 


6565 


45 


0.56 


H/3C 


4540 


114 


0.53 
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Table 4. Total exposure times in seconds per filter for NGC 4321 and NGC 4254. The aperture of the telescope used is given 
in brackets 



Filter NGC 4321 NGC 4254 



Hq 


7200 (0.8 


m) 


1800 (2.5 


m) 


RaC 


5400 (0.8 


m) 


1800 (2.5 


m) 


H/3 


5400 (1.5 


m) 


1800 (2.5 


m) 


H/3C 


3600 (1.5 


m) 


1800 (2.5 


m) 


B 


7200 (0.8 


m) 


1200 (2.5 


m) 


I 


7200 (0.8 


m) 


600 (2.5 


m) 
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Table 5. Limiting fluxes and magnitudes for NGC 4321 and NGC 4254, at three times the r.m.s. of the sky background 



Filter 


NGC 4321 


NGC 4254 


Ha 


1.1 • 10"" erg cm"^ s"^ 


3.4 • 10"" erg cm~^ s'^ 


H/3 


5.1 • 10-1* erg cm-2 s'^ 


6.0 • 10"" erg cm'^ s'^ 


B 


23.34 mag arcsec"'^ 


25.40 mag arcsec"^ 


I 


22.18 mag arcsec"^ 


23.10 mag arcsec"^ 
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